Scalable Ellipsoidal Classification for Biparite Quantum States 
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Abstract 

Tie Separability Problem is approached from the 
perspective of Ellipsoidal Classification. A Density- 
Operator of dimension N can be represented as a vec- 
tor in a real vector space of dimension N 2 — 1, whose 
components are the projections of the matrix onto 
some selected basis. We suggest a method to test sep- 
arability, based on successive optimization programs. 
First, we Gnd the Minimum Volume Covering Ellip- 
soid that encloses a particular set of properly vector- 
ized bipartite separable states, and then we compute 
the Euclidean distance of an arbitrary vectorized bi- 
partite Density Operator to this ellipsoid. If the vec- 
torized Density Operator falls inside the ellipsoid, it 
is regarded as separable, otherwise it will be taken 
as entangled. Our method is scalable and can be 
implemented straightforwardly in any desired dimen- 
sion. Moreover, we show that it allows for detection 
of Bound Entangled States. 

1 Introduction 

Entanglement is an essentially new resource for com- 
munication that lies at the very heart of Quantum 
Information Theory. Among several other exploits, 
it has been shown to permit information transfer at 
higher rates than classically expected, thank to the 
existence of quantum correlations between separated 
parties [1] |2j E]. Thus, a theory of Entanglement 
which offers a qualitative description as well as quan- 
titative measures is highly desirable. 

For pure states, a satisfactory theory exists and 



there are procedures both to detect and quantify the 
entanglement of a given state. For instance, the 
Von Neumann Entropy of any of the reduced Den- 
sity Operators gives a measure of the entanglement 
contained in the whole system, being zero if and only 
if the state is separable, and is directly related to 
the Schmidt Decomposition of the Density Operator. 
For mixed states of dimension less or equal than 6, 
i.e. for 2 x 2 and 2x3 systems, Positivity under Par- 
tial Transposition (PPT) is a necessary and sufficient 
condition for separability [HIS]. However, in the gen- 
eral case, things are more involved, and it has been 
shown that, for dimension greater than 6, the prob- 
lem of checking whether a mixed state is entangled 
or separable is NP-hard [6]. Thus, one should not 
expect to find an exact solution to the Separability 
Problem. 

In the recent years some approximations based on 
different relaxations of the problem have been sug- 
gested. Doherty et aZ.[7J|H] proposed a test based on 
nested semidefinite programs to search for symmet- 
ric extensions of a given bipartite system, which is 
a stronger condition for separability than the PPT 
criterion. Brandao and Vianna showed that the Sep- 
arability Problem can be cast as a robust semidef- 
inite program, which is an NP-hard problem, and 
proposed both probabilistic and deterministic relax- 
ations [TU]. Note that these tests only provide 
bounds to the certainty of the answer to the ques- 
tion "is the state p separable or entangled?" , but they 
don't offer a quantification of the entanglement con- 
tained in p. 

Our method relies on the concept of the Mini- 
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mum Volume Covering Ellipsoid (also known as the 
Lowner-John Ellipsoid) of a given set of vectors in 
vector space R N _1 , isomorphic to the space of Den- 
sity Operators of dimension N. We show that, for 
vectors representing separable bipartite Density Op- 
erators, it is possible to identify this ellipsoid with 
the boundary of the set of bipartite separable states, 
up to some accuracy. Interestingly, this approach can 
be used to estimate the amount of entanglement con- 
tained in p. 

In next section, we briefly review fundamental is- 
sues concerning geometric entanglement measures. 
In section 3 we introduce our method and apply it 
to 2 x 2 and 2x3 systems. In section 4, we show that 
it allows for detection of Bound Entanglement. 

2 Review of Geometric Entan- 
glement Measures 

The study of entanglement from the geometrical 
point of view offers insightful ways both to detect and 
to quantify quantum correlations between different 
systems. First one is known as Entanglement Wit- 
nesses (EW), which split the space of states in two 
parts in such a way that the set of separable states 
lies strictly inside one of the resulting half-spaces. 
We denote the set of all Density Operators by T>, 
and the set of separable states by S. It follows from 
probability theory that both sets are convex [TT] . Let 
W = be an EW, then: 

(W, p) < 0, for some P eV\S (1) 

(W, a) > 0, Vcr S S (2) 

Clearly this defines a hyperplane dividing T>: 
(W,a) > (W,p). Here (.) is the inner product for 
matrices (A,B) = Tr{A^B). This product naturally 
defines the Frobenius norm \\A\\p = ^Tr(A^A) 

Another measure is the Hilbert- Schmidt Distance, 
Dhs, of a Density Operator to the set of separable 
states S: 

D HS = min \\p - a\\ F (3) 



Note that computing Dhs involves the obtention 
of cr = Ps(p), the projection of p onto S. An Opti- 
mal Entanglement Witness, that is, an EW for which 
([I]) holds for the largest number of states in T> \ S, 
will be tangent to S at this point. 

2.1 Projection of Density Operators 
onto a Vector Space 

Any Density Operator p of dimension N can be writ- 
ten in vectorized form: 

j N 2 -l 
i=l 

where 1 at is the unit matrix of dimension N and 
{(Tj}^ _1 is some basis of self-adjoint, trace-free ma- 
trices, such that: 

Trto-iO-j) = aSij (4) 
N 

n = —Tri^p) (5) 
a 

and 7 is chosen to normalize r. The Frobenius 
norm induces an Euclidean norm in R N -\ For pure 
states, ||r|| 2 — 1, while for mixed states ||r|| 2 < 1. 
This illustrates the convexity of the set of states. 

Thus, a one-to-one relation between Density Oper- 
ators of dimension N and vectors in R w _1 can been 
established. This allows us to switch from the matrix 
space picture to the real vector space picture, which 
will prove useful for our method. 

2.2 Duality between Hilbert-Schmidt 
Distance and Entanglement Wit- 
nesses 

It is a general result from geometric optimization that 
the problem of finding a separating hyperplane be- 
tween a point p and a convex set C is dual to the 
problem of finding the distance between C and p [T^] . 
In matrix space language, this duality can be illus- 
trated as follows. Let p and S be an arbitrary Density 
Operator and the set of separable states. The prob- 
lem ^ can be posed as: 



2 



An ellipsoid centered at x c can be expressed as: 



min |M|f (6) 
such that p — a = t 
a e S 

with variables r and a. The Lagrangian of ^ is: 

C=\\t\\f + (W, P -<t-t) (7) 

where W is the Lagrange multiplicator associated 
to the equality constraint. It is not hard to see that 
it represents a a hyperplane. Noting that (W, r) < 
II^II-fIMI-Fj the dual function can be written as: 

g(W)= mm [(W,p)-{W 1 a) + \\r\\ F (l-\\W\\ F + S)} 

(8) 

where the parameter 5 > is related to the relative 
orientation between the hyperplane represented by 
W, and the line going from the separable set S to the 
Density Operator p, and it is equal to zero if and only 
if they are perpendicular. For ^ to be bounded from 
below in r we must include the additional constraint 
- S < 1 . So the dual problem of ^ is: 

max . [xom[(W,p)-(W,a)]] (9) 

||VC]| — 5<1 c£5 

It is straightforward to check that the optimal 
value of (J9j) is attained if and only if W is an Opti- 
mal EW. This result, also known as the Bcrtlmann- 
Narnhoffer-Thirring Theorem (see Ref. [13]), will let 
us trace a link between the entanglement detection 
and quantification problems (compare also with Refs. 

MM) 

3 Ellipsoidal Classification 

The basic premise of our method is that the set of sep- 
arable states S can be approximated by a Minimum 
Volume Covering Ellipsoid (MVCE) of an ensemble 
of vectors corresponding to some separable Density 
Operators. Then, the following classification scheme 
can be adopted: if a vector falls inside the MVCE, it 
will be taken as separable, and if it falls outside, it 
will be regarded as entangled. 



£ = {x|(x-x c ) T A(x-x c ) < 1} (10) 

where A — A T is a positive definite matrix of di- 
mension N 2 — 1 . The volume of this ellipsoid is pro- 
portional to det(A -1 / 2 ). 

Since in matrix space quadratic forms are not 
defined, we need to work in a real vector space 
to build this ellipsoid. We first obtain an en- 
semble of "separable vectors" by means of tenso- 
rially multiplying states along all directions spec- 
ified by some canonical basis. For instance, in 
the 2x2 case, this ensemble would be {x| ep } = 
{(1,0,0) <g> (1,0,0), (1,0,0) <g> (-1,0,0), (1,0,0) <g> 
(0,1,0), (1,0, 0) (8(0,-1,0),. ..,(0,0, -1) (8(0,0,-1)} 
. Later on we will see that it is convenient to vary 
the norm ||x^ ep ||2 of these vectors. This procedure 
ensures that all vectors will lie as spaced as possible 
in the separable set S. Secondly, we minimize the 
volume of an ellipsoid, constrained to have all gener- 
ated "separable vectors" falling inside it. One way to 
obtain the MVCE of this ensemble would be to solve 
the following problem: 

min logdetA -1 (11) 

such that (xf p - x c ) T A(xf p - x c ) < 1 

with variables A and x c . Here, logarithm was taken 
in order to drop off proportionality terms. Despite 
the exponential growth of the dimension of the as- 
sociated vector space, interior point methods used 
for minimization still converge polynomially to a so- 
lution in dimension as large as 1000, or more [T2"] . 
Throughout this work, we used the Disciplined Con- 
vex Programming Rule-Set implemented in the Mat- 
lab Package CVX [16]. 

3.1 Results for 2x2 and 2x3 Systems 

The Separability Problem is solved for 2 x 2 and 2x3 
systems, thanks to the PPT criterion. We use this 
fact to benchmark our method. For a sample of 1000 
"separable vectors" and 1000 "entangled vectors" , we 
solved two different problems: first one is the original 
problem ([3} casted as: 
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Figure 1: The vertices of the polytope are the gener- 
ated "separable vectors" of which we find the MVCE. 
The larger set corresponds to the whole space of Den- 
sity Operators 





2x2 Syste 
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False Negatives 


0.1 


962 
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868 
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687 
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484 
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287 


15 
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180 


184 
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92 


410 


0.8 


32 


600 


0.9 
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755 


1.0 





873 



Table 1: Number of misclassified vectors in a sam- 
ple of 1000 "separable vectors" and 1000 "entangled 
vectors" , as a function of the Euclidean norm of the 
vectors of the separable ensemble 



mm ||p — u\\f 
such that o Ta > 



(12) 



which gives the true results. The variable is a, and 
u Ta denotes Partial Transpose in subsystem A. The 



second problem was to find the MVCE through (11) 



and compute the distance to this ellipsoid in a similar 
way: 

min 1 1 r — s 1 1 2 (13) 

such that (s — x c ) T ^4(s — x c ) < 1 

where r and s stand for the vectorized counterparts 
of p and a. The results obtained for pure vectors 
(||x| ep || 2 = 1) are rather discouraging: whereas none 
of the generated "separable vectors" fell outside the 
MVCE, only 12.7% of the "entangled vectors" were 
detected. This led us to shrink our ellipsoid by re- 
ducing the norm of the generated ensemble {x^ ep }. 
At the expense of letting some "separable vectors" 
fall outside the ellipsoid, we increase the number of 
correctly classified "entangled vectors" . This recalls 



a binary hypothesis test, and we shall adopt the cor- 
responding terminology to expose our results. The 
event that a true "separable vector" falls outside the 
MVCE will be a false positive, while if an "entan- 
gled vector" falls inside the MVCE, it will be false 
negative. Stepwise reducing the norm of the vectors 
belonging to the separable ensemble we obtained Ta- 
bles 1 and 2. 

There is a trade-off between the number of cor- 
rectly classified states and non-ambiguousness of the 
test. The relevant area of 2 x 2 systems is between 
norms ||x| ep || 2 = 0.6 and ||x* ep || 2 = 0.5, as can be 
seen in Fig. [2] A measure of entanglement should be 
as unambiguous as possible, and thus we argue that 
the best choice is ||x| ep || 2 = 0.5, since for this case 
only about 1.5% of the "entangled vectors" are mis- 
classified. For this choice, in general, a vector will be 
misclassified 15.1% of the time. For 2x3 systems (see 
Fig. [3]), the MCVE approximates somewhat less effi- 
ciently the separable set. However, still 76.8% of the 
vectors are correctly classified, in the area comprised 
between ||xf p || 2 = 0.5 and ||xf p || 2 = 0.4. In these 
systems, it misclassifies at least 5.2% of the "separa- 
ble vectors" . 
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FALSE NEGATIVES 





2x3 Syste 


ms 


Norm 


False Positives 


False Negatives 


0.1 


949 
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812 
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597 
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427 


52 
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269 


196 
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160 


389 
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80 


572 
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34 


699 
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11 


807 


1.0 





900 



Table 2: Number of misclassified vectors in a sam- 
ple of 1000 "separable vectors" and 1000 "entangled 
vectors" , as a function of the Euclidean norm of the 
vectors of the separable ensemble 



Figure 3: False Negatives versus False Positives for 
2x3 systems. The error probability can be reduced 
to 23.2%, between ||xf p || 2 = 0.5 and ||xf p || 2 = 0.4 



3.2 Pseudo-Entanglement Witnesses 
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For a vector space endowed with the Euclidean norm, 
there is a simple way to construct a tangent hyper- 
plane to a given ellipsoid. We can use this fact to 
build realistic observables amenable to a laboratory 
setting. 

The tangent hyperplane to the ellipsoid can be ex- 
pressed as: 



V x [(x - x c ) T A(x - x c ) - l] SD (r - s ) = (14) 

where Sq = Ps(r) is the projection of the vectorized 
Density Operator under study onto the MVCE. It can 
be expressed in affine form as: 



Figure 2: False Negatives versus False Positives for 
2x2 systems, showing that there exists an area where 
the probability of wrongly classifying a vector can be 
brought down to 15.1%, between ||x* ep || 2 = 0.6 and 



0.5 



(s - x c ) T A(r -x c ) = 1 



(15) 



(compare with Ref. |17j). It is important to keep 
in mind that, although the hyperplanes introduced in 



(15 1 very much resemble an Entanglement Witness, 
they are not so in general. This is because our MCVE 
may in general be a proper subset of the separable 
set S, and no tangent hyperplane to this MVCE will 
strictly separate S from any entangled state. Nev- 
ertheless, these Pseudo-EW can be used to estimate 
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10 20 30 40 50 60 70 80 90 100 

100 RANDOMLY SELECTED "ENTANGLED VECTORS" 



Figure 4: For 100 random "entangled vectors" of a 
2x2 system, the continuous black line is the true 
distance to the separable set S, while the dashed line 
stands for the distance of the vectors for a MVCE 
of "separable vectors" of norm ||x^ ep ||2 = 0.5. At 
the bottom, the pointed line represents the distances 



obtained for norm llx 



sep I 



= 1 



the amount of entanglement contained in a given en- 
tangled matrix p via ([9]), which at the optimal value 
will be equal to (|6]l|12]~ For an illustration of entan- 
glement estimation sec Fig. |4j 



4 Bound Entanglement Detec- 
tion 

For composite systems of dimension higher than 6, 
there is a special kind of entangled states that can- 
not be used, in principle, to enhance communication. 
The entanglement contained in these states cannot 
be distilled to obtain pure entangled states [T5] , and 
it receives the name of Bound Entanglement (BE). 
The PPT criterion fails to detect this kind of states, 
and it just becomes a necessary condition for quan- 
tum correlations to arise. Other criteria to detect 
and quantify BE have been proposed, such as non- 
decomposable EW [TO], Schmidt number Witnesses 
|20j . and, more recently, a geometric approach based 
on separating hyperplanes [21] . 
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x 3 Systems 


Norm 


Detected States 


0.1 


1000 


0.2 


1000 
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1000 


0.4 


1000 
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1000 
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1000 
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226 


0.8 


149 


0.9 


107 


1.0 


79 



Table 3: 1000 bound entangled states were generated, 
with parameter "a" running from 0.001 to 1. The 
distances of the associated vectors to the different 
MVCEs were obtained. For norms of the separable 
ensemble of 0.6 and below, all bound entangled states 
were detected 

Our approach is in the spirit of the latter of the 
aforementioned methods, but instead of hyperplanes, 
we shall use the MVCE in order to detect BE. In- 
tuitively, the ellipsoid covering a set of "separable 
vectors" should leave bound entangled states on its 
outside. This fact is studied in 3X3 systems, where a 
parametrization of bound entangled states, due to P. 
Horodecki, is available [32]. These states pbe depend 
on a scalar a £ [0, 1], and are given by: 



Pbe (a) 
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Surprisingly, for norms of the generated separable 
ensemble of 0.6 and below, all bound entangled states 
are detected. The obtained results are shown in Table 
3. 

As expected, the distance to the MVCE of the de- 
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Figure 5: There is a linear dependence between the 
distance to the MVCE and the norm of the associated 
Density Operator 



tected states linearly depends on the norm of the as- 
sociated Density Operator. This interdependence is 
depicted in Fig. [5] 

Although, for sufficiently small norms, all BE 
states are detected, one should keep in mind that, 
in these cases, an unacceptable number of false pos- 
itives is introduced. This fact precludes a truly re- 
liable criterion. On top of this, not all types of BE 
states were tested, so no complete evidence for de- 
tection of BE is provided. Nevertheless, this gives a 
quantitative picture to the intuitive idea that, since 
BE states lie outside the separable set, some of them 
will necessarily be detected by our MVCE method. 



5 Conclusions and Further Per- 
spectives 

Our method presents some drawbacks. First impor- 
tant one is ambiguousness, since false positives must 
be necessarily introduced, to some extent. Also, the 
efficiency of the MVCE decreases as 1/JV as the di- 
mensionality of the problem grows. This can be com- 
pensated reducing the norm of the "separable vec- 
tors" belonging to the generated ensemble. Although 



the reduction of the norm is geometrically justified, 
we are still lacking a qualitative physical argumenta- 
tion, if there is any. 

On the other hand, it presents at least two appeal- 
ing advantages. Interior point methods behave well 
in very high dimensions, so in principle any discrete 
quantum system can be accounted for. For this rea- 
son the method is scalable. Moreover, relatively few 
resources are needed in contrast to other methods, 
and in fact, once the MVCE is obtained, a decision 
is made in very few seconds. 

This method could probably be refined allow- 
ing some robust quadratic classification, or using a 
polyhedral approximation instead of computing the 
MVCE. However, this has the taste of being NP-hard. 

The author acknowledges fruitful conversations 
with B. Beferull, F. Brandao and J. de Vicente. 
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